Feedback Loops Between Fields and Underlying Space Curvature: an Augmented 

Lagrangian Approach 
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We demonstrate a systematic implementation of coupling between a scalar field and the geometry 
of the space (curve, surface, etc.) which carries the field. This naturally gives rise to a feedback 
mechanism between the field and the geometry. We develop a systematic model for the feedback 
in a general form, inspired by a specific implementation in the context of molecular dynamics (the 
so-called Rahman-Parrinello molecular dynamics, or RP-MD). We use a generalized Lagrangian that 
allows for the coupling of the space's metric tensor (the first fundamental form) to the scalar field, 
and add terms motivated by RP-MD. We present two implementations of the scheme: one in which 
the metric is only time-dependent [which gives rise to ordinary differential equation (ODE) for its 
temporal evolution], and one with spatio-temporal dependence [wherein the metric's evolution is 
governed by a partial differential equation (PDE)]. Numerical results are reported for the (1-f 1)- 
dimensional model with a nonlinearity of the sine-Gordon type. 
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Recently, much attention has been focused on soft- 
condensed-matter objects, such as vesicles, microtubules, 
and membranes [1-4]. Many nanoscale physical systems, 
including nanotubes and electronic and photonic waveg- 
uide structures [5,6], have nontrivial geometry and are 
influenced by substrate effects. These classes of systems, 
many of which are inherently nonlinear, raise the ques- 
tion of the interplay between nonlinearity and a substrate 
with variable curvature. Of particular interest is a pos- 
sibility of developing curvature in the substrate due to 
forces generated by the nonlinear field. The resulting 
curvature can in turn affect the field. 

There is an increasing body of literature dealing with 
the interplay of nonlinearity and a curved substrate. Usu- 
ally, however, the substrate geometry is assumed to be 
fixed, see, e.g., [7]. Nevertheless, for many applications, 
ranging from condensed matter to optics to biophysics, 
it is relevant to introduce models that admit a flexible 
substrate, which is affected by the field(s) that it carries, 
as well as feeding back into the field dynamics. In this 
situation, equations for the fields in a nonlinear system 
abutting on the flexible substrate should include both the 
field dynamics proper and the feedback coupling to the 
substrate. Equations for the evolution of the substrate 
should in turn be affected by the evolution of the field. A 
prototypical physical example of this type is Euler buck- 
ling [8], where the evolution of a thermal profile causes 
the underlying surface to buckle (and hence locally mod- 
ify its curvature). 

In a discrete setting, a model of this type has recently 
been presented in [9]. However, it was limited to a sys- 
tem of masses coupled by nonlinear springs. Some studies 
have also been performed in a special case of the contin- 
uum limit of classical spin systems (such as the Heisen- 
berg chain) coupled to the curvature; geometric frustra- 



tion was found to arise in such settings [10]. 

About twenty years ago, a problem similar to the 
theme of our study was examined in the context of molec- 
ular dynamics studies of structural transitions in crystals. 
In particular, in a series of papers [11], Rahman and Par- 
rinello introduced a new idea for studying such transi- 
tions by means of an augmented Lagrangian that would 
account for the degrees of freedom of the "box" (the cell) 
in which the MD particles lie. In studying the time evolu- 
tion of the box dynamics (naturally obtained through the 
Euler-Lagrange equations of the augmented Lagrangian 
for the box degrees of freedom) , they were able to identify 
structural transitions (under external shear) from square 
to hexagonal patterns, fee to bee etc. We will hereafter 
refer to this technology as the RP-MD method. The rele- 
vant Lagrangian for the particles and the box in this case 
reads: 

L^\Y. m,s,Gs, - V{r,,) + \wTv{f^ f), (1) 

where is the mass of the i-th particle, is its vectorial 
velocity, the spatial part G of the spatiotemporal metric 
tensor may be represented in terms of another matrix / 
as G — f"^ f {G is positive definite). and Tr denote 
the transposition and trace respectively, rij and V are 
distances between the particles and the potential of in- 
teraction between them, and W is an effective mass of 
the box. 

Our purpose in this work is to extend the RP-MD 
methodology to the case of a continuum scalar field, cou- 
pled to either a spatially averaged geometric characteris- 
tic ("average curvature"), which will give rise to an ODE, 
or to a spatiotemporal curvature field, that will generate 
a PDE. The continuum field may represent, e.g., a chem- 
ical concentration propagating over a membrane, or a 
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salt solution, causing the swelling of a polymer gel [12]. 
or an envelope wave of the electric field in nanosystems. 
The spatiotemporal metric is assumed to have the simple 
form, 
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Wc will first consider the general case, where G is a d x d 
matrix, d being the space dimension. 

One can define a field-type generalization of the RP- 
MD model, with a scalar field 6, as follows: 
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where V{(f)) is a potential governing the nonlinear evolu- 
tion of the field </> and / = f{t) (only). If / (and hence 
G) is a function of both spatial coordinates and time, 
an clastic-energy term [13] should be added to the La- 
grangian (3), so that it becomes 
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An objective of this brief report is to propose two mod- 
els that include the feedback to the curvature, and, si- 
multaneously, admit as particular solutions the (unper- 
turbed) solutions for the flat (original) metric. In partic- 
ular, wc choose G = 1 + ./^ (so that the metric is positive 
definite and for f — has as a special case the original 
Minkowskian metric). Notice a slight deviation in our 
choice from the RP case of G = [14]. Although the 
formulation of Eqs. (3)-(4) is very general, we hereafter 
focus on the d = 1 case that we will examine in more 
detail. 

Assuming initially that / = f{t) only (e.g., including 
only an effect of the "mean curvature" on the scalar-field 
dynamics), the Lagrangian (3) becomes 
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(5) 



Then, the resulting equations of motion (to which we will 
hereafter refer as model A) are 



<l>tt = (1 + - 

Wftt = -f j c^ldx, 



dV 
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(6) 
(7) 



where the subscripts stand for the corresponding partial 
derivatives. 



Notice that the fimction / is directly related to the 
scalar curvature of the 1-d space. In particular, the Ricci 
scalar, which is i? = 2R\'2i2 1 det{g) [16] in the general 
case, in the 1-d case is i? = —Ijajj, where / = yl + Z^. 

On the other hand, for a metric with both spatial and 
temporal dependence (e.g., for / = f{x,t)), one arrives 
at the following Lagrangian: 
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The ensuing coupled equations for the scalar field and 
the curvature (to which we will refer as model B) are 



</>«= ((l+/')0x),- 

Wfu = wu, - f<f>l. 
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(9) 
(10) 



As a particular application of models A and B, we ex- 
amine the physically ubiquitous sine-Gordon (sG) poten- 
tial, V{(t)) = 1 - cos((/)) [15]. It is clear that Eqs. (6)-(7) 
and (9)-(10) have particular solutions with / = 0, for 
which the latter equation of each pair is satisfied triv- 
ially, while the former reduces exactly to the sG equation. 
Basic solitary-wave solutions of the sG equation are the 
topological soliton (kink). 



i>ii{x, t) = 4 tan ^ [exp (7(0; — xq — vt))] 
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where v is its velocity, 
Lorentz factor, and xq 
the kink's center, and 



4 tan 



7 = (1 - t,2)-i/2 is the 
is the initial position of 
the breather, (/)br(a;, t) = 
u'^)/u}'^ 8uv[u^{t — v{x — Xq))] X 



scch [7-\/l — uj'^{x — xq — vt)\ }, where uj is the frequency 
of its internal oscillations (0 < a; < 1). 

The results of the interaction of the kink with the cur- 
vature in model A are shown in Fig. 1 [17] . The curvature 
variable /, initialized with a small random value, per- 
forms smooth oscillations with a frequency of cj w 2.866. 
Notice that this is natural in this case, since the kink 
has an approximately fixed "mass", Mk = u^dx, 
which can be found to be Mk = 8.247 for the velocity 
V « 0.25). Then, Eq. (7) predicts the frequency of these 
oscillations ui « V^k =2.871, which is very close to the 
above-mentioned numerically exact value. Fast small- 
amplitude oscillations of the kink's velocity, observed in 
Fig. 1 both with and without the curvature, are due to 
"hopping" over sites of a lattice (with spacing h = 0.1) 
employed in the numerical scheme which solves Eq. (7). 
Notice, however, that in the top panel the mean veloc- 
ity is w 0.2503, while in the bottom panel it is w 0.2499, 
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hence the curvature oscillations increase the kink's veloc- 
ity. This may be anticipated due to the presence of the 
positive definite factor 1 + in front of (j)xx in Eq. (6), 
which is expected to renormalize . 

The curvature-breather interaction in model A is 
shown in Fig. 2. The frequency of the breather does 
not change significantly (it fluctuates between 0.87 and 
0.93), but its amplitude decreases substantially (by more 
than 50%), resulting in its becoming much more mobile 
(the velocity increases to w 0.425 from the initial value 
0.25). 

In model B, we examine collision of a kink with a lo- 
calized pulse of the substrate field /. For the breather, 
we have obtained results which are qualitatively similar 
to those presented below for the kink. We create the cur- 
vature pulse to the left or to the right of the kink. As (j)^ 
vanishes far from the kink, the equation (10) for / be- 
comes a linear wave equation. Hence, we observe splitting 
of the pulse into left- and one right-traveling ones. In the 
case where the kink is initially to the left of the pulse, it 
collides with the left-propagating fragment of the (split) 
pulse. In the opposite case, the kink is eventually caught 
by the co-propagating right fragment of the (split) pulse. 
Numerical simulations shown in Fig. 3 demonstrate that 
the collision with the counter-propagating pulse reduces 
the kink's velocity, while the interaction with the co- 
propagating pulse gives rise to an increase of the velocity. 
In particular, in the former case, the mean speed of the 
kink after the collision is « 0.2476, while in the latter 
one, it increases to « 0.2527. Notice that in both cases 
a small fragment of the pulse that collides with the kink 
passes through it, while a larger fraction of the pulse is 
reflected by it. The latter feature may be explained by a 
momentum-balance analysis. 

In conclusion, we have extended field theory in the 
spirit of the Rahman-Parrinello Molecular Dynamics 
technique. The resulting equations couple the spatio- 
temporal evolution of the field to that of the underlying 
curvature of the space which carries the field. Coupled 
equations for the temporal or spatiotemporal evolution of 
the metric are obtained in a general setting, and, as an 
example, are solved together with the sine-Gordon field 
equation. The purely temporal evolution of the metric 
has been found to increase the velocity of the field soli- 
tons, while the model allowing spatio-temporal evolution 
of the metric can induce both increase and decrease of 
the velocity. 

It would be particularly interesting to extend the mod- 
els A and B to higher dimensions. It is also worth study- 
ing how the local evolution of the curvature affects kine- 
matics and dynamics of the solitons, and to correlate such 
observations with the behavior of reactant chemical con- 
centrations in chemical or biological environments with 
non-trivial geometry [18]. 
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FIG. 1. The top left panel shows the kink's spatial profile at t = 60 in model A. The initial kink profile (centered around the 
new center position) is shown by the dashed line, and is practically indistinguishable from the solution at i = 60, indicating 
that the kink maintains its shape. The bottom left panel shows smooth oscillations of f{t). The top and bottom right panels 
show, respectively, the kink's velocity v{t) vs. t, and the same quantity but for / = 0. 
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FIG. 2. The top left and right panels show, respectively, the breather in model A at the end of the simulation period, t = 120, 
and its position vs. time (solid), as compared to that which it would have moving at the initial velocity, v = 0.25 (dashed). 
The middle left panel shows the time evolution of f{t), while the middle right and bottom left panels present exchange between 
the curvature-mode's energy, Ef = + if /A) J 4>ldx, and the rest of the energy, = E — Ef. Finally, the bottom right 
panel shows the time evolution of the field 4> at the center of the breather. 
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FIG. 3. The left part of the figure shows spatial profiles of the fields <j) / at t = 60 in model B. It also shows the time 
evolution of the energies — E — Ef and Ef — J {l/2){ft + + (f4>xY /2)dx, as well as the kink's velocity as a function of 
time. The kink's center is initially at a; = 90, while the center of the curvature pulse is at a: = 110. The right subplots show 
counterparts of those in the left part, but for the case of the kink and curvature pulse initially centered at a; = 110 and x = 90, 
respectively. 
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